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Abstract 

The lift on a strongly non-spherical vesicle in a bounded shear flow, is studied in the case the 
membrane moves with a velocity, which is a linear function of the coordinates. The magnitude of 
the induced drift is calculated as a function of the axes lengths, of the distance from the wall, and 
of the ratio of the cell to the solvent viscosity. It appears that the main mechanism for lift, in the 
presence of tank-treading motions, is the fixed orientation of the vesicle with respect to the flow. 
Tank-treading vesicles in suspensions, flowing through narrow gaps under small Reynolds number 
conditions: Re = Lv/v < 1 (with L the gap width, v the flow velocity and v the viscosity) migrate 
away from the walls, with a velocity that is 0{Re~^) larger than predicted by inertia. 
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An important property of suspensions, both in Couette and channel flow conditions, is the 
ability to change their effective viscosity, by concentration of the particles, near the centre of the 
flow 1^. One of the most striking examples is the behavior of red cells in small blood vessels, 
where this phenomenon takes the name of Fahraeus-Lindqwist effect Q . 

In the case of spherical particles, it is well known that the only mechanism, which is able 
to cause a transverse migration of particles, is the presence of inertial corrections to the Stokes 
equation, which describes the response of the solvent to the particles ^, ^. This is not true 
anymore, however, when spherical symmetry is lost. In this case, a particle can behave like a sort 
of wing and move in a direction different from that of the flow. 

Very often, in particular when one deals, as in the case of red cells Q, with vesicles, or otherwise, 
with droplets of an unmiscible liquid, the particle shape is approximately that of an ellipsoid. If the 
ellipsoid is rigid, it will not maintain a fixed orientation, but will rotate as a whole in an unsteady 
fashion, carrying on a kind of flipping motion One suspects in this case, that spherical symmetry 
is in some way recovered and that Stokes dynamics is not able to produce by itself a transverse 
drift. Things change in the case of an object with a compliant structure, where tank-treading 
motions become possible and a fixed orientation ensues. 

Recently, the behavior of a spheroidal vesicle, undergoing a tank-treading motion in a bounded 
shear flow, has been analysed using perturbation theory around the spherical particle case . For 
a spherically symmetric geometry, a Green function for the Stokes equation is available and this 
equation can then be solved, for arbitrary boundary conditions on the surface of the cell. 

For strongly non-sphcrical objects, however, the only tractable situation is that of the rigid 
ellipsoid, studied by Jeffery [||. Keller and Skalak |]lO| showed that Jeffery's analysis could be 
generalised to situations, in which the boundary condition on the cell surface, is that the velocity 
be a generic linear function of the coordinates. The special case they considered, with a cell 
assumed to have an ellipsoidal shape, and with the points on its surface moving with constant 
angular velocity, provided a useful approximation for the locally area preserving motion of a red 
cell membrane. 

Of course, this approach solves only a part of the problem, leaving out the determination of 
the shape and internal motions of a cell in an external shear. A self-consistent description, like the 
one contained e.g. in jl^, 13|, is however not really necessary here, at least for large distances 



from the walls, in which case, no large deformations are expected. 

In this paper, the theory of Jeffery, and of Keller and Skalak, is used to calculate the lift on a 
vesicle, undergoing the special kind of tank-treading motion studied in . Analysis carried on in 
M, showed in the perturbative regime, that different kinds of tank-treading motions produce lifts, 



which are essentially of the same order of magnitude. (More precisely, the case considered in |10 
gave a lift, which was approximately twice as large as that of a cell, with a membrane moving in 
an exactly area preserving fashion). 

As in , the analysis is confined to a range of distances from the wall, such that a calculation of 
the cell- wall interaction, based on the method of images is possible, and that the far field expression 
for the velocity disturbance from the cell can be used. Extrapolating the results to the region near 
the wall, it appears that a strongly non-spherical tank-treading vesicle, migrates towards the centre 
of the fluid, with a speed of the same order of the difference of velocity, due to the external shear, 
that exists across its body. This speed is typically much larger than the one produced by inertia. 

Consider a cell, whose shape is that of an ellipsoid, with axes 02 < ai < as, immersed in a shear 
flow, which, in an appropiate reference system {xi, X2, X3}, translating with the ellipsoid and with 
origin at its centre, can be written in the form: v = KX2e3. Indicate with {x[, x'2, x'^} the reference 
system in which the equation for the ellipsoid surface can be written in the diagonal form: 

/ 2 / 2 / 2 
\ + \ + (1) 
at ai 

In a state of tank-treading motion, the cell will keep a fixed orientation, with x[ = xi, and x'^ 
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forming an angle 6 with respect to x^, with < ^ < 7r/4 for k > (see Fig. 1). For vanishing 
internal viscosity, or in the almost spherical case, 9 —> 7r/4, corresponding to the strain part of the 
shear flow dominating the dynamics and forcing the ellipsoid and the strain axes to get aligned. 
In the other limit, 9 — signals the crossover to the vorticity dominated, flipping motion regime. 




Figure 1: Orientation of a tank-treading cell in a plane shear flow v = KX2e3. 

The constant angular velocity tank-treading motion, that is taken into exam, is described by 
the equation for the velocity at the membrane [[L0[ : 



02 
03 



/ / , as 

02 



2'=3 



(2) 



This kind of motion can easily be shown not to be area-preserving by an amount of the order of 
log 03/02; however, Eqn. (2) goes beyond providing an expression for which allows analytical 
calculations. Recently, Kraus et Al. have obtained semi-analytical expressions, for 9 and oj in 
terms of vb{xi = 0), and for vb{xi — 0) itself, keeping into account all the couplings between the 
membrane and the fluid. It turns out that, in the spheroidal limit, the expression given by Eqn. 
(2) for vb{xi = 0), (and then also for 9 and u) is equal to the one obtained in Jilt . 

The boundary condition for the velocity perturbation v produced by the cell, is given by 
V = V — vs on the cell surface, and can be written in the form: 



0, 



V', = K{bx', + if- 04), 4 = + if + 04), 



(3) 



where: 



sin 261 



and f = - 
■' 2 



are the components of the strain matrix, while 



cos2t^ + L0[ 

03 02 



(4) 



(0.2 

Vos 



03 
02 



(5) 
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is the vorticity of v. Jeffery's theory |g] gives the following far field expression for the velocity 
disturbance: 

v^-^^<i>(x) (6) 
where $ = Ax'^ + Bx'^ + Cx'^ + 2Fx'2x'^ and R — (010203)^ . In the unprimed laboratory frame: 

^ = Ax\ + {^^-^ cos2e + - F fiiTi2e^xl 

^ " ^. cos 29 + + F sin 2^1 xi 



2 2 
+ [(B-C)sin20 + 2Fcos2^]a;2a;3 (7) 
The coefficients A, B, C and F are given by the expressions [|j: 

^ = z ' ^ = z ' ^ = z = 

^=^5-T74-^ ^here: Z = eK'a^ + a^'a^' + a^), (») 

where a; = R~^ai, and the a-coefficients, which depend solely on the geometry of the cell, are 
given by: 

, _ f"° dA , dA 

"i^X MTa 

and: 



(ai + A)(a2 + A)A' ^ y„ + A)(a? + A)A ' 



where 



A = ^{al + X){al + X){al + X). (10) 
Writing ai = a2 = (1 — e)^i?, 03 = (1 + e)5_R and taking the limit e ^ 0, one obtains: 

$ ~ sin 2^[-2ix? + (i - F)a;^ + (i + F)a;^] + 2{F + F cos 29)x2X3, (11) 



F = -^ A = -^, and ^.-(l-^). (12) 



where: 

5e ^ 5e , 5 5e ' 

16 ^ "112' 

Upon substitution into Eqn. (6), Eqns. (11-12) give the expression for the velocity disturbance by 
a spheroidal particle, obtained in 

The orientation angle and the tank-treading frequency u) were calculated in ||l^ by Keller 
and Skalak, imposing angular momentum conservation, and requiring that viscous dissipation be 
equal to the work done by v on the system: 



cos 20 = 



1 



(13) 
(14) 



zi^-l — I, Z2 = a^{a^+a2), 23 = \ • (15) 

2 V a2 a3 / a2 03 



2 (l-e^)z2-2 

_ cos2^ 

~zi((l-e^)z2-2)' 

In the equations above: 

1 /as a2\ , ,.2 , »2n 03 , 02 
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and is the ratio of the inner to outer viscosity, with the inner viscosity being in general a weighed 
average of the viscosity of the membrane and that of the fluid inside pO| . The situation in which 
the magnitude of the RHS of Eqn. (13) is greater than one (as it happens for instance in the rigid 
particle limit — > oo) corresponds to the flipping motion regime. 

The derivation of Eqns. (13-15) becomes very simple when — 1 and the membrane viscosity 
is equal to zero, which corresponds to the case considered in pl| ]. It is necessary to know the 
expression of the velocity inside the cell, but it is easy to show that this is given by the same Eqn. 
(2) describing the membrane motion (see e.g. |^ and the expression for the inner solutions to the 
Stokes equation derived there). Imposing that no energy be dissipated in the membrane, forces the 
vorticity to be continuous across the cell surface, leading to the condition: ^ = 0, while, requiring 
that no torque be present in Eqn. (3), leads to the strain being aligned with the ellipsoid axes: 
/ = 0. These two conditions imply: V2 = Kbx'2 and Wg = —nhx'^, and using Eqns. (4-5), the result 
of Eqns. (13-15) for = 1 is recovered. 

A plane wall, placed perpendicular to the a;2-axis, at a distance / from the cell, causes a 
disturbance v' in the velocity field, which, at the cell centre results in a net drift away from 
the wall. This drift can be calculated, if the wall is sufficiently far away, by the image method, 
imposing the no-slip boundary condition at the wall: = —v. The calculation, whose details are 
illustrated in |9|, is carried on introducing a representation in terms of potentials: v' = V0-f V x A, 
V • A = 0, and Fourier transforming in xi and x^. The potentials satisfying the Stokes equation 
and incomprcssibility for v' are: 

Ak(a;2) = (2:2 - 0Akexp(fc(a;2 - 0), (/)k(a;2) = exp(fc(x2 - /)). (16) 

where fu{x2) — J dxidx^ exp(— ik • x)/(x), / = {ip, A} and k = {fci, k^}. Imposing the boundary 
condition on vj^(a:2), by means of Eqns. (6-7), and inverse Fourier transforming the result, the 
large distance expression for the drift = —^2(0) is obtained: 

v^-^ (17) 

where U = U{tv, ^1,^2), with = Ui/a^, is twice the coefficient of the X2 term entering $, as given 
in Eqn. (7): 

[/ = (B - C)cos2^+ (S + C) - 2Fsin2^. (18) 

In the almost spherical limit, Eqns. (11-12) lead to the result of ~ i|esin2^, with ^ ~ 7r/4. 
The quantity U, which is of order 1 for large departures from sphericity, is a dimensionless velocity, 
normalised to kR, and gives the magnitude of the drift close to the wall. 

Given Eqns. (6-9) and (13-15), the functional form of U{e„,ri,r2) can easily be explored. As 
it was to be expected, there is a strong correspondence between U{tu,ri,r2) and 9{ev,ri,r2)] in 
particular: U{6 0) = 0, which corresponds to the lift vanishing at the transition to flipping 
motion. 

If both parameters and r2 are kept fixed, the strongest lift is achieved when ri = r2 (see 
Fig. 2), which parallels the effect, observed in ||l^, that tank-treading motions are more robust in 
prolate rather than oblate ellipsoids. This is not due, however, to the orientation angle 9, being 
the largest in one of the two cases. It appears instead that is independent of ri: — 9{e^, r2). 

The strongest variations in U takes place, when either the viscosity ratio e^, or the non-sphericity 
(which can be parametrized by r^^) are large. This corresponds to the intuition that strongly non- 
spherical, stiff cells, prefer to stay in a flipping motion state. In Fig. 3, as an example of strongly 
non spherical object, the red cell model studied in pO| is considered, and compared with the case 
of a spheroid. The red cell, both in a stretched and an unstretched state, make a transition to 
flipping motion roughly at = 3, while the lift on a spheroid, remains almost constant up to 
en ~ 10. For fixed e^, the lift grows with the degree of non-sphericity, from [/ = at e = 0, and 
reaching a maximum at a critical e = ^{tv), when ri = r2. For e > e, U decreases and reaches zero 
at the transition to flipping motion (see Fig. 4). 
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Figure 2: Plot of the drift U in function of the eccentricity e = (y^^)^, for = 1, in the case of 
axisymmetric elHpsoids: (a) prolate ri = r2; (6) oblate ri = 1. The orientation angle 9 (c) remains 
the same in the two cases. 
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Figure 3: Dependence of U on the viscosity ratio e^; red cell at rest: (a) ri = 1, r2 = 0.286; (6) 
stretched red cell: ri = r2 = 0.213; (c) a prolate ellipsoid: ri = r2 = 0.5; an oblate spheroid: 
n = 1; (d) r2 = 0.8. 
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Figure 4: Dependence of U on e, in the case of prolate ellipsoids, for different values of e,y. Notice 
the maxima moving towards e = as the viscosity ratio is increased. 



For e^, small, it appears from Fig. 2, that the variations in U/e, that take place in the case of 
axisymmetric ellipsoids, are of the same order of those due to changes in the type of tank-treading 
motion. Thus, for small e^, the linear theory in keeps doing an acceptable job also for large 
values of e. Notice however, that the linear theory remains unable to make predictions on 6, which 
instead, is kept fixed at ^ = tt/4. 

In the case of the red cell model studied in jl^ , Fig. 5 illustrates the dependence of the lift on 
the degree of stretching and on e^. It must be remembered that, in that case, the ratio e,^ must 
account for the viscosity of the membrane as well as that of the hemoglobin inside . However, 
for e,y ~ 1, corresponding to red cells in a strongly viscous solvent, one sees that U ~ 0.3, with not 
much variation when the cell is stretched by the shear. One can compare the prediction of Eqn. 
(17), with the contribution to lift from inertia, which is: foi^^R ^ with v the viscosity of the 

suspension, and /o of the order of a few tenths Q. Thus, if one has a suspension of red cells in a 
viscous solvent, forced to flow in a gap of thickness L, the condition for the lift on the red cells, to 
be dominated by inertia, becomes kL^ > i.e. a large Reynolds number condition for the flow in 
the gap. 

From an experimental point of view, this means that the time T ^ L/v^ for cells to concentrate 
near the centre of a channel flow, as revealed e.g. through viscometric measurement, will obey 
different laws in the two regimes, so that T = T{k) must have a kink at k ~ i/L~^, with T oc k for 
K < R and T a for k > R. These experiments are complicated however by the contribution of 
vesicle deformability to the suspension effective viscosity [|l5| . 

The goal of this research was to extend the analysis carried on in to a non-perturbative 
regime. It appears that the drift produced by tank-treading motions approaches, close to the wall, 
the magnitude of the velocity difference at the particle scale kR. It must be stressed, that in this 
region, the approximations leading to Eqn. (17) break up, so that the expression for the drift 
this equation provides, has a value only as an order of magnitude estimate. For larger values of 
l/R, however, Eqn. (17) becomes exact. This furnishes then, a test for numerical solvers, once a 
boundary condition like that provided by Eqn. (2) is imposed. For ? ~ i?, a more quantitative 
analysis would require a full solution of the fluid equations, in the presence of both the cell and 
the wall. Due to the complicated geometry of the problem, a numerical approach would become 



8 



0.5 
0.45 

0.4 
0.35 

0.3 

U 0.25 
0.2 
0.15 
0.1 
0.05 


0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

r2/ri 

Figure 5: Values of the lift on a red cell in solvents of different viscosity, in function of the stretching: 
r2/ri. 

necessary. 

Looking at the problem from a different angle, analysis carried on in suggested a rather 
strong dependence of the drift on the exact nature of the tank-treading motion, and similarly, on 
deviations from a strictly ellipsoidal cell shape. The magnitude of the drift appears to be strongly 
influenced, also by the viscosity of the membrane and that of the fluid inside the cell. Hence, 
considering both measurement difficulties, and how details in shape and membrane motion are 
affected by perturbations, an exact knowledge of the value of U close to the wall is perhaps not 
very interesting, and the use of a simple expression for v^, like the one provided by Eqn. (2), is 
appropiate. 

As an aside, these results imply that, a cell capable of producing tank-treading motions, has at 
its disposal a propulsion mechanism, which, beyond obtaining most of the necessary energy from 
the external flow, can be controlled by small changes in cell shape. A similar mechanism for cell 
propulsion has been suggested in 

Despite these strong dependencies, one expects the production of lift from the fixed orientation 
and non-sphericity of the cell, to be in itself a robust phenomenon. In particular, the same effect 
should be present, as assumed earlier, even when the large distance condition is relaxed, and more 
important, when quadratic corrections in the shear are taken into account. (In fact, in Poiseuille 
and channel flow conflgurations, one would expect a stabilising effect on the cell orientation 
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